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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
er eH NOT BE CONSTRUED AS A COMMIT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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++ 
FACILITY: MATH LIBRARY 
ABSTRACT: 
This module contains the routine MTHSAMOD: 


It returns the remainder of the division of argl/arg2 using 
the following equation: 


arg] - (int(argl/arg2))*arg2 


; AUTHOR: Bob Hanek, CREATION DATE: 21-DEC-1982 


; MODIFIED BY: 
Jeffrey C. Wiener, 29-DEC-82 
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-SBTTL HISTORY ; Detailed Current Edit History 
Edit History for Version 3.0: 
3-001 Original version of a complete rewrite JCW 29-DEC-82 
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ht a DECLARATIONS BSEp= 1984 11:90:52 EMTRTL. SRE IATHAMOD MAR: 1 

«-SBTTL DECLARATIONS 
INCLUDE FILES: 

NONE 
EXTERNAL SYMBOLS: 
-DSABL GBL 

-EXTRN MTHSSSIGNAL 

EXTRN ATHSKTINVARGMAT 
LIBRARY MACROS CALLS: 

SSFDEF ; Define SFS (stack frame) symbols 
EQUATED SYMBOLS: 

NONE 
OWN STORAGE: 

NONE 
PSECT DECLARATIONS: 

-PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 
CONSTANTS: 

NONE 


; Force all external symbols to be declared 
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0 ~SBTTL MTHSAMOD - F REAL*4 remainder 
: FUNCTIONAL DESCRIPTION: 


wr 


FUNCTION VALUE: 


Remainder of the division of argl/arg2 is returned as an 
F_floating point value. 


IMPLICIT OUTPUTS: 
NONE 
COMPLETION CODES: 


94 : Return the remainder of argh args in Part point format 
a6 3 Remainder = argi = (int(argl/arg2))*arg 
$3 ; The algorithm used to evaluate the AMOD function is as follows: | 
6 : X = the first argument. 
1 8 Y = the second argument. 
1 § step 1. m = the exponent of Y. 
1 § 3 n = the exponent of X. 
80 104 : ite <0 end with result = X | 
0000 105; step 2. I = the fractional part of xwg024 
4 196 3 J = the fractional part of Y¥*2*24 
00 107 ; If 1 >= Teled 
385 198 3 step 3. c =c- 
00 109 ; If c_<_0 go to step 7. 
000 110; step 4. L = 2*31* 
4 41 3 Il=el- 5 anche) 
3 c=ce 
0000 118 3 If ¢ >= ° go to step 4. 
0000 114; step 5.c=c + 31 
0000 115 ; If c_>= 0 go to step 7. 
44 116 3 step 6. L = 2*c#l 
000 117 ; I=L = Jeint(L/J) 
00 ; 118 3 step 7. Result = 2*(m-24) * | 
0 0 0 ; CALLING SEQUENCE: 
8 0 ' ¢ : Remainder.wf.v = MTHSAMOD (dividend.rf.r,. divisor.rf.r) 
0 4 > INPUT PARAMETERS: 
0 6 1 8 : The two input parameters are F_floating-point values. They are 
0 4 : 3 passed by reference. 
00000004 eis. % : DIVIDEND = 4 ; Dividend = X in the algorithm. 
00000008 0 : ? DIVISOR = 8 ; Divisor = Y in the algorithm, 
: § IMPLICIT INPUTS: 
134 NONE 
: 5 
Hi 
if 
i 
168 
144 
145 
146 
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NONE | 
SIDE EFFECTS: 
MTHS_INVARGMAT = Invalid argument to math Library if the divisor is zero. 


MTHS" FLOUNDMAT - eles ing underflow in math library is signaled if 
the FU bit is set in the callers PSL. 


-ENTRY MTHSAMOD, “M<CR2, R3, R4> 
52 08 BC 50 MOVF +3445) aha R2 : Rg = Y, the divisor 
52 8000 8F AA BiCw2 #*x8000, R2 3 R2 = {Y¥} 
10 BNEQ START 7 ivi <0 
50 01 or ASHL #15, #1, RO : t¥t=0. Division by zero case 
7E O0O°BF 9A MOVZBL #MTHSK_INVARGMAT, -(SP) t error code 
00000000'GF 01 CALLS #1, G*ATHSSSIGNAL ; signal the error 
50 04 BC 00008000 8F cB START: BICL3 #*XxX8000, @DIVIDEND(AP), RO ; RO = {KX} 


54 52 FFFFOO7F 8F CB BICL @*XFFFFOO7F, R2, R4 ; R4 = m the exponent of Y 
53 50 eid ed P. BICL #*XFFFFOO7F, RO, R3 ; R3 = nm the exponent of Xx 
5 ae ° SUBL2 R4, R3 3; R3 = c¢ = men 
; plus some low order bits 
69 «#19 BLSS GET_SIGN : If c<O then {Xi > {Y¥! and the 
$ result is X 


OOOO CCCOOCCoOoCoooooo 


STEP_2 
Extract the fraction part of X*2°24, called I, and the 
fractional part of Y¥*2*24, called J. 


WAAAAANI 
ow 


point representation, the hidden bit needs to be added Tnto the 
internal representation since the number is to be converted to 
an integer value. 


So 

oO 

onmon— _ 

ow o 
SOSOSCOSOOOSSOSOSSOSOSOOOOSOSOOOSOOOOOOOOOOOOOOOOOOOO 


Oooooooooooeoo 


FDO DIDI et st rs ee ed nd od od od od 

SOOCOOOOOOOO OO O09 09 09 09 09 09 09 09 NINN NSIS NS NSN PAA AAA AAA A NTU 

WN OS OONA UNE WIN 9 OD NAME WIN 0 ODNOU EW)  O OO WO UNE A OO OO NOUS WH" OWO@on~n 
= 


B 
38 
50 7F80 8F AA 38 BICW #*Xx7F80, RO ; Clear the exponent field 
50 90000080 , a 40 ADDL2 #*x80, RO : Replace hidden bit 
5 0 10 rh ROTL #16, RO, RO ; Convert to integer (RO = I) 
52 ree FF OAA 4 BICW #*X7F80, R2 ; Clear the exponent field 
52 900099 0 8F co 30 ADDL2 #*x80, R > Replace hidden bit 
5 :..) C of ROTL #16, R2, R2 ; Convert to integer (R2 = J) 
52 50 01 22 CMPL RO Re 3; Compare | and J 
03 19 BLSS STEP : Branch if I < J 
50 2 C2 ? SUBL2 R2, RO gs 5 g€ee fe J 
6 3+ 
3 : STEP_3 
6 3 Convert c = exponent of X = exponent of Y into an integer. 


After the exponent bits are removed from the internal F_floating 
| 
| 
| 
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MTHSAMOD = F REAL*4 remainder 
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| 
Subtract 31 from 5, in order to Goternine if an iteration | 


of the algorithm needed. If c-51>=0 then go to STEP_S. 


ROTL #-7, R3, R3 : Convert § to an integer vglue 
SUBW #31, R : Check shift count, c = c=31 
BLSS STEP_5 ; Branch, if ¢ < 0 


STEP_4 
Compute 1 = L = Jeint(2*c*#Il/J) by rem(2*c*l, J) since I and 
J were scaled to integer values. 


ROTL 


#-1, RO, R1 RO/R1=2°31*1. This and the next 
BICL3 #*X?FFFFFFF, R1, RO 


two gna trys Shans are equivalent 
to ASHQ #31,R0,R0, but are faster 
RO/R1 contains L = 2*31#] 

RO = rem(2*31*I,J) 

Check shift count, c = c-31 
Branch if c >=0 


Sete Se Se Ge Ge Se 


E Ro. RO R1, RO 
SUBW2 #4  Rigy 
BGEQ STEP_4 


#31, R3 Restore shift count, c = c+#31 
STEP_7 :; If zero, branch to STEP_7 


STEP_6 
Compute I = L = Jeint(2*c*I/J) by rem(2*c*l, J) since I and 
J were scaled to integer values. 


LRL R1 : 

SHQ R3, RO, RO ; RO = 2*c#l 

DIV R2, RO, R1, RO 3; RO = rem(2*c*l, J) 

VILF RO, RO : Convert I to floating point 
ete, #*x4C00, RO ; RO = me a * | 

DDW 4, R ; RO = 2*(m-24) * | 

LSS UNDERF LOW ; Branch if underflow occured 
STW @DiVIDEND (AP) 3; Check for sign of result 

GE RETURN and adjust answer accordingly 


Q ; 
EGF RO, RO 3 


RO ; set up default result to 0.0 
BBC #SFSV_FU, SFSW_SAVE_PSW(FP), NO_FU 
; Branch if caller has not enabled F 


PUSHL #MTHSK_FLOUNDMAT ; Report MTH$_FLOUNDMAT 
CALLS #1, G*ATHSSSIGNAL : si nal the érror 
: Return 


IVIDEND = "F 
IVISOR = 8 
GET E¥ AG R 
MTHS$SIGNAL eeeceere §6X 
MTHS RG 


eeeeeere 8 § X 
eeeeeeee Xx 


Initialization 
Command processing 


Pass 
Symbol table sort 
Pass 

Symbol table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


The working set Limit was 1050 


Page faults 


5219 bytes (11 pages) of virtua 

There were 10 pages of symbol table space allocated t 

source Lines were read in Pass 1, produc ing 13 object records in Pass 
pages of virtual memory were used to define 7 mac 
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Allocation PSECT No. Attributes 

00000000 ( 0.) 00 ¢ OQ.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
00000000 0.) O1¢ +1.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
000000C ( 194.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
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ap 


ages. 
onery were used to buffer the intermediate code. 
hold 44 non-local age 0 local symbols. 


ros. 
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MTHSAMOD 
VAX=11 Macro Run Statistics 
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! Macro Library statistics ! 


tere rt ences mene oe seoen ocean accent 


Macro Library name Macros defined 
~$255$DUA28: (SYSLIBISTARLET.MLB;2 eager 
88 GETS were required to define 4 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHAMOD/OBJ=OBJ$:MTHAMOD MSRC$:MTHAMOD/UPDATE=(ENH$:MTHAMOD) 
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